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If the choice of generalized coordinates for determining 
the transient response of a non-symmetric structure were not 
eigenvectors but were modes of deformations due to operating 
loads, there would be certain advantages. Among these would be: 
1. the economy of requiring only a small number of modes, 2. the 
avoidance of having either to cull out certain non-participating 
modes or to retain the non-participating modes at the expense of 
having to operate with larger order matrices, and 3. the con- 
fidence of getting well converged solutions. Using load response 
modes as generalized coordinates is properly classified as the 
Ritz Method. 

The interest in the case of non-symmetric stiffness 
derives from structures with active control systems. The assym- 
metry comes from the sensor being situated at a different loca- 
tion than the actuator. Loads that would be typical of those 
used in the design of control systems are: externally applied 
forces and pressures, vernier jets whose firing is commanded by 
the control system, and constraints from appendages assigned to 
attach points. Such a solution method is developed here as a 
DMAP ALTER packet to the Statics Rigid Format in NASTRAN. The 
Inertia Relief Rigid Format would appear to be a more natural 
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host for such an ALTER packet, because it carries out full proc- 
essing of mass. Its prohibition against boundary constraints 
disqualifies it, hence Statics with modifications to process mass 
is used as host. 

THEORETICAL APPROACH 

The theory will be organized in 5 parts. It will de- 
scribe the mathematical decisions on using the original stiffness 
matrices in the development of fundamentals first. It will 
describe the generation of harmonics from each fundamental, 
second. Then it will develop the adjoint vectors third. Having 
a full complement of primary and adjoint generalized vectors, 
the next step is to orthogonalize them and finally to integrate 
them into an actual solution by constructing the generalized mass 
and generalized stiffness and reconciling the form of the gener- 
alized damping. 

The methods presented herein are an outgrowth of a new 
non-col located sensor actuator analysis method under development 
by H. P. Frisch at the NASA/GSFC, Code 712.^ The motivation for 
the NASTRAN/DMAP implementation presented herein is to provide a 
working capability which can be used for both current practical 
applications and for the evaluation of the new analysis method 
prior to its inclusion into the general purpose multi-flexible 
body data preparation program FEMDA. 


1. Frish, H.P. "IAC Program FEMDA, Theory and User's Guide, 
Interface from Structural Analysis Output Data to Input Data for 
Multi Flexible Body Dynamics Analysis," NASA Tech Brief Draft, 
June 13 , 1988 . (Call author for status info 301-286-8730) 
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FUNDAMENTAL MODES 


A structure must be defined according to its elastic 
distribution, its damping distribution, its mass distribution, 
its boundary conditions and its complement of those loading 
conditions which are active during its operation under a control 
system. Each such loading that can be applied independently of 
some of the other loadings should be treated as a distinct defor- 
mation-producing condition; i.e. as one producing a unique static 
response shape. Within the theory of linear elasticity the 
magnitudes of static response to a static loading varies directly 
as the magnitude of the loading, so a simple unit magnitude of 
load is all that is necessary to establish the shape of an in- 
dividual mode into which the structure deforms. A sketch of a 
cantilevered beam with an end load illustrates the point. The 
loads are graduated from 1 unit, to 2 units, to 3 units. 



The ratio of deformation to load stays constant at 
L^/(3EI> as the load varies. The same holds true at positions 
other than the tip deflections. For instance, the deflection at 
an interior point, such as L/s where 1 < s < °°, is 
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1 ' 

— J , therefore 

j-j , which shows that the ratio 

of deformation to load at an interior point stays constant as 
load varies. The notion of shape that is independent of ampli- 
tude can be easily depicted by a sketch of a violin string being 
played in its fundamental mode. 




All three show the string sounding the same pitch ( frequency ) but 
at different loudnesses (amplitude). All 3 deformations are 
considered to have the same shape. 

Therefore, in order to get started, a controlled struc- 
ture must be exercised with a unit load for each loading condi- 
tion, which can vary in magnitude independently of other loads. 
This will produce, what will be called, the set of fundamental 
modes for the complement of loading conditions. It was assumed 
at the start. of this derivation that the set of loads produced a 
set of unique static modes. This assumption needs to be tested 
against a criterion for uniqueness. The criterion that is ger- 
main here is linear independence. The modes need to be put on an 
equal footing for such a test by normalizing them uniformly. 
There is a choice of methods for normalizing at this point. 
Because this is a dynamics problem, one would be inclined toward 
mass normalization, but for now a simple Euclidean length will 
suffice to put the modes on an equal footing to establish linear 
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independence. Later, when the modes are orthogonal i zed, they 
will be normalized to mass. Normalizing to maximum or to the 
amplitude of a common position will be excluded, because of the 
bias that would be introduced in the linear independence check. 

Initially static modes are to be obtained by solving n 
loadings . 

(1) [ K Ll] [{ u 1 } ' { u a} { u n}]t, = [{ P l}’{ P 2}' { P n}] L 

■■■> [“uJKJ-hJ- 

Each of the ju^j of CU3 will be individually normalized by its 
Euclidean length. Compute the individual normalizing constants 
according to 

<2, [{«/{«,}] - n r 

Each term of |u^J will now be divided by the square root of 
the normalizing constant n^. Name this normalized fundamental 

(3) PHI sub j, i.e. ~zz{ u -j} = 

In order to use these static deformations as generalized 

coordinates, they should be linearly independent. Now we are 

faced with the decision as to what criterion of linear independ- 
ence to use and what tolerance to allow. The classical defini- 
tion of linear independence is that the Gramian > 0. The 
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Gramian is a determinant of the matrix of the dot products of 
each coordinate into every other coordinate. Translated to the 
context in which we want to consider it, assume that the group of 
vectors 
them i s 


♦l ♦„ 


( 4 ) G( <(>^ <(> 2 r 
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If G ( <t>^ <f> n > = 0, the set of <|) n are a dependent set, but if it 

is > 0 the set of 4» n are linearly independent. In a finite 
number system, such as that under which a digital computer oper- 
ates, the establishment of a true zero is difficult, if not 
impossible. Consequently the decision as to what criterion to 
use rests on the practical consideration of how much impurity to 
allow in the set that we want to use as the expansion functions. 
If the Gramian is just slightly greater than zero, it implies 
that yes a functional relationship can be set up for one in the 
set with respect to the others, because imperfections creep in to 
contaminate the zero computation. Now if all <j>'s were normalized 
to unity and constituted a truly linearly independent set, then 
the value of the Gramian would be unity. 

G( unity normalized (J>) = 1. 


2. TENSOR ANALYSIS, by I. S. Sokolnikoff; pp 6; John Wiley, 1951 
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Any departure from 1 in the Gramian of this unity normalized set 
means that some imperfection is creeping in and the further the 
descent from 1 towards zero implies more and more imperfections 
away from true linear independence. So What! Our goal is to 
have a good set of vectors so that when we expand our solution in 
them, we will get good accuracy and get good convergence. If we 
have a set of linearly independent vectors but too few of them to 
span the range of actions in which our structure will operate, we 
will fail to converge close enough to a correct solution. If we 
admit too many vectors that are almost independent, but do have 
some imperfections, the answers will contain biased emphases and 
will definitely contribute amplitudes that are too great in some 
modes. In further consideration of the practical factors that 
will govern our decision, we ask, “What will it cost to find out 
if there is linear independence?" In the case of the Gramian, 
the computations involve the creation of a full matrix of vector 
dot products that must then be decomposed and finally the product 
of all diagonal terms of the decomposed matrix forms the Gramian 
determinant. Decompositions are expensive, so a method other 
than using the Gramian would be worth while to investigate. 

Another approach is to look at the ingredients of the 
Gramian i.e. the individual matrix dot products. By definition 
the dot product of 2 vectors A and B is 

A . B = | A | x | B | cos 0, where 0 is the generalized angle 
between the 2 vectors. This can be extended to vectors in N-di- 
mensional space. Set up a criterion based upon the size of the 
angle that a trial vector <{> makes with each of those (j> n that have 
already been judged linearly independent. 
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(5) 6. = cos 

k,n 

threshhold value. Ideally x would be ir/2 for an orthogonal set 

with perfect linear independence. The poorest possible value 

would be that for which <J> and <t>. are coincident, i.e. zero 

n k 

angle. This criterion can be rephrased by saying that it tests 
how well cos 0 compares with cos ir/2 
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t, where x is some 
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7 ) 


The desire is to hold the angle between test vectors to be some- 
where between a threshhold and tt/ 2. If the threshhold angle were 
tt/3, k = cos tt/3 = .5, then the test would require the cosine of 
the angle between test vectors to be less than 0.5 : 



0.5 = k for all 1< n < k. 


Once k has been decided upon, the test is carried out against 
every i(i n tor n < k. 

In the case of the cosine test, the matrix of dot pro- 
ducts would have to be formed as in the case of the Gramian, but 
no decomposition need be done. The absolute value of each term 
is compared to k. 

The Gramian test has the advantage of making its decision 
by comparing only one number against a threshold while the cosine 
test involves comparison of every ratio in a column against a 
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threshhold. The term by term processing involves only simple 
operations; the net result is that the cosine test is much less 
expensive than the Gramian test. The cosine test has been chosen 
for this method. 

har mon i c s 

Modes, determined solely on the basis of static loads, 
are questionable to apply without supplement to the solution of 
dynamics, because they are devoid of inertia effects. Supplemen- 
tal modal vectors can be generated by finding the deformation due 
to forces derived from accelerating masses distributed through 
the structure by an amplitude equal to the vector of elastic 
deformation. Call the deformation from inertia effects, ac- 
celerated by amplitudes derived from the fundamental mode, the 
first harmonic. A second harmonic can be generated by the scheme 
used to generate the first harmonic, except that the the inertias 
are now accelerated through amplitudes derived from the first 
harmonic. Similarly, a third harmonic can be generated from the 
acceleration of mass through the amplitudes of the second har- 
monic, etc. Eventually the upper harmonics will tend towards 
congruence, so there will be an nth harmonic beyond which no 
distinct modes will be added. The measure to be used for finding 
the useful limit will be linear independence. 

The set of linearly independent modes, consisting of a 
group of fundamentals plus groups of harmonics associated with 
each fundamental, when normalized to the Euclidean length, then 
orthogonal ized, will be used as modes of generalized coordinates 
in expanding the behavior of a structure under the management of 
a control system. The mathematics of these inertia modes fol- 
lows. The mass of the structure |~M r T 1 will be accelerated by an 
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amplitude distributed spatially according to the shape of the 
normalized - but not orthogonalized - fundamental mode j- to 
create an inertia forcing for the first harmonic |F^^|; i.e. 


( 7 ) 

( 8 ) 


Ki} • KJk) - uhere 

/ • • % , d^ \ ' 2 \ 

\^f/ = {^f — “ cos = cos wt )| . 

d t 


The shape of the deformation through which the mass will be 

accelerated is established by <|>~. The effect of the term 

2 1 
(uj cos wt) is to merely amplify the shape as a function of time. 

Our interest, at this stage of the derivation, is only in the 

shape and not the total dynamic response, therefore the forcing 

can be treated as a statics problem with the spatial distribution 

of the set of accelerations, limited to any instant of time, 

therefore 

I > » | | 

<|>£ a 4)j and the resulting static force is 


i9 > Ki} - KJk}- 


Apply this force to the structure and solve for the response. 


<10) [ K Ll] { u hl } { F hl} [ M Ll] {^f } ' 

from which the response can be explicitly isolated: 


,n > ki} - kikkiJk}- 

The so obtained will be normalized by the Euclidean length 
and will be tested for linear independence , which if accepted, 
will be named j^i} and augment the complement of Ritz 
modes. This first harmonic will have a shape of deformation 
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sufficiently unique to make it worth while to consider it as the 
basis for accelerating the mass through its spatial behavior to 
obtain a second harmonic, similar to the way that the fundamental 
was used in the generation of the first harmonic. Mathematically 
the method of forming the second harmonic follows the pattern 
already established for the first harmonic. Form the forcing 


( 12 ) 


Ka} ■ Kl]{U- 

2 

Let {<J> hl | = {♦hi} 


(cos wt) and limit the value to its 


static amplitude | <t> hl | - | 4» hl | / then solve for the static 

response ju^J from 

,13 ' [ K Ll] ( u h2> - K 2 } - KlJKi} • 

Extract | u h2 } explicitly. 

<14) { U h2} = [ K Ll] [^Ll] {^ hl} * 

The |u^ 2 | so obtained will be normalized by its Euclidean 
length and will be tested for linear independence, which if ac- 
cepted, will be named an< * augment the complement of 

Ritz modes. A question arises as to the extent to which a can- 
didate harmonic should be tested for linear independence. Should 
it be tested against every other vector established up to this 
point, or should the candidate harmonic be tested for linear 
independence only against its parent fundamental and the 
harmonics that are spawned from that fundamental alone? From an 
algorithmic standpoint the latter route is favored, because all 
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quantities stay within an inner loop. The merit of this deci- 
sion can be tested and be replaced if need be. It is known that 
as the recursion steps are carried out for higher harmonics, the 
deformations will tend toward congruence so there is a definite 
need for testing each new candidate harmonic against its parent 
and siblings. If there is no physical risk for so limiting the 
linear independence check to the family associated with just one 
fundamental, it will be opted for here. 


A certain pattern starts to appear from the development 
of these two harmonics. The matrix product 1 [ M Ll] is use<i 
repeatedly. Consequently, it can be generated once and saved for 


recall in the generation of any level of harmonics of any fun- 
damental. Name this matrix product 


(15) 


CS0Li3 



where i can take on the BCD 


character for the primary or the adjoint mode, and * can take on 
either blank for primary or T (for transpose; for the adjoint. 
Discussion of adjoint mode generation will be taken up subse- 
quently. 


Capitalizing on the pattern that has been revealed, all 
first harmonics for all fundamentals can be generated in a series 
of matrix operation as follows, as adapted for the primaries: 
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(ISA. [{u 1 } > {u 2 } i {u 3 } i .... i {u 1 } r ... i {u n }] sh = 

[V] 3 ' [ M Ll] [{*l} ,{* 2 } ,{* 3 } K} . ' ' ' . (♦»}] Sf 

This can be compressed into 
( 15B ) [U H1 ] s = [K L J s [m ll ] [» f ] g 

Equations <15,15A,15B) represent all possible inertia response 
raw data for forming first harmonic modes of primary fundamen- 
tals . 


Normalization of these responses can also be performed on 
all vectors treated as a matrix. The Euclidean length can be 
extracted from the multiplication of hi] by itself . 

T 

(iso [u H1 ] [u H1 ] - [u H1 sg], 

Strip off the diagonal, take the inverse of each, fol- 
lowed by its square root, to form a diagonal matrix of scale 
factors mode by mode. 



n I - * 5 


1 15D) 

L' U H1 S0 v] — 

L scal hi] 


Apply this matrix of scale factors to the first harmonic respons- 
es, of fundamental inertia loadings, as a post multiplication 
operation to get a set of candidate normalized harmonic modes. 

USE. [u H1 j [SCAL H1 ] = [♦“]. 
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This matrix of candidate first harmonics needs to be tested for 
linear independence. Testing will proceed in two parts. All 
first harmonics will be tested against all other first harmonics, 
then all first harmonics will be tested against all fundamental 
modes. Start with the harmonics by themselves. The cosine test 
involves talcing the dot product of every mode against every other 
mode. This is done in a single matrix operation. 

««*> [*Kf[C] ■ Ku] • 

Examining the jjxiT^jJ matrix closely, one can recognize that the 
first row represents the dot product of the harmonic (first in 
sequence) against each of those in the set of harmonics. The 
second row represents the dot product of the harmonic (second in 
sequence) against each of the set of harmonics. Et cetera. 
Consequently, the next step is to strip off one row to examine 
how well this candidate harmonic holds up in the cosine test 
versus other first harmonics. Next another row is stripped off 
and this candidate is tested and so on until all candidates have 
been examined. If the vectors had not been initially normalized 
to length, it would have been necessary to do so at this point to 
form the cosines. As a consequence, the matrix DOT^ consists of 
all cosine terms. Going back now to the first row, some detail 
will reveal a pattern for systematizing all of the candidates. 

Select one term at a time starting with the 2nd and take 
its absolute value then compare that value with k. Shift the 
index to the 3rd and do the same. Continue until either a value 
greater than k is encountered or until the the end of the row is 
reached. If any term tests greater than k, catalog the row 
number and proceed to the next row. All successful candidates 
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will be eligible to be tested against the matrix of fundamental 
modes. The harmonic vs. fundamental test will patterned after 
the harmonic vs. harmonic test. All successful candidate 
harmonics will be held in reserve to form the basis of 2nd 
harmonics before they will be merged with, but in sequence after, 
the fundamental modes. 

The number of successful first harmonic modes may be 
fewer than those in the set of fundamentals. This does not 
matter, because the method of computing harmonics is independent 
of the size of the order of the vectors from which they are 
derived. The generation of second harmonic modes and higher will 
proceed along the pattern just outlined for first harmonic modes, 
except that after the responses to inertia loads have been com- 
puted, the successful modes from which they were derived, will be 
merged into the matrix of previous Ritz modes. Now the i— 
set of harmonics will be tested for linear independence vs. not 
only themselves but against all other Ritz vectors including 
fundamentals and all previous order successful harmonics up 
through the (i-1) — . 

The generation of higher harmonics will be subject to a 
choice of two limitations. The analyst may want to limit the 
maximum number of harmonics to be admitted for any particular 
investigation because of, say, a study in a low frequency domain. 
He can invoke such control by giving the value of the maximum 
number of harmonics to the parameter MODSPEC. Declining to 
assign a value to MODSPEC will cause the number of harmonics to 
be limited by those that pass the linear independence check. 
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ADJOINT MODES 


For the non-symmetric problem, adjoint vectors are re- 
quired to obtain the reduced order coefficient matrices. There 
are no set rules for introducing the adjoint basis. In the 
spirit of Lanczos method, a trial method is introduced and 
refined. For lack of any better trial scheme, the starting 
matrices of the adjoint system that will be used here will 


LL 


consist of the transpose of the original 
[ K Ll] ' (whictl wi H also be non-symmetric), and the original 
vectors 


K 


J matrix; i.e. 

load 


The static solution of this adjoint system under the n 
original loads yields a set of responses | v j| : 

<16) Kl][KM v 2 } { V n}j ' [{ P l}'{ P 2} { P n}] 

----> fr&JfcJ ■ kJ • 

Each of the of CV3 will be individually normalized by its 

Euclidean length. Compute the individual normalizing constants 
according to 


,17) [{ v j) W] = V 

Each term of j v j| will now be divided by the square root of the 
normalizing constant a ^ . Name this normalized fundamental 

< 18 > PS I sub j, i.e. -~(vj s { 41 J. 
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Each adjoint fundamental will be checked for linear independence 
against all other adjoint fundamentals. If any adjoint fundamen- 
tal mode fails the linear independence check, and if its primary 
companion passed, both primary and adjoint will be discarded. 
This is necessary in order to retain a uniform sequence when 
setting up generalized mass and generalized stiffness. The 
cosine test will be used to certify linear independence. Iner- 
tial harmonics of adjoint fundamentals will be generated in the 
same manner as those of the primary fundamentals. Once again the 
linear independence of the adjoint harmonics will be checked only 
against its parent and siblings, instead of the currently estab- 
lished set of ALL Ritz modes. 

ORTHOGONALIZATION 


The solution of the differential equations is enhanced if 
the generalized coordinates used to span the response space are 
orthogonal. Our set of linearly independent vectors can be 
orthogonalized. At this point we have a pair of bases vectors 
£03 and £0} that are each separately linearly independent. After 
orthogonalizing they will be given the symbols ££3 and £fl} re- 
spectively. Options for orthogonalization and for associated 
constraint weightings were given due consideration. Before 
deciding on what options to choose, it will help to review the 
ultimate application. 

The dynamic equation in metric coordinates is: 

• • » 

(19) CMJtx 3 + CBJCx} + CK3Cxj = F(t). 
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Transform the metric coordinates to the primary bases 
vectors as generalized coordinates. Assume that the orthogonali- 
zation of <J> into Q has already taken place. Let 

(20) £x J = C£Dt£3 , giving 

CMIICCKS 3 + CB3CCKU + tKDCUtU = F(t). 

At this point the columns only of the coefficient matri- 
ces M, B, and K have been transformed to generalized form. Next 
the rows of the matrices must be transformed. In the symmetric 

3 

case this is done by using the transform of the symmetric 
modes. However, in this the unsymmetric case, we pre-multiply by 
the transpose of the adjoint bases vectors. Assume that the 
orthogonalization of into il has taken place. 

(22) Cii3 T CM3CC3£i’i + Cii3 T CB3CC3£U + O^’cKJCCHU = Cii3 T F(t). 

The desire is to decouple the equations as much as possible. To 
start with, we want the generalized mass to be a unity matrix; 
i.e. 

(23) Cil3 T CM3CC3 = 1 \] - 

Constraining the generalized mass to unity will affect the nor- 
malization to such an extent that the generalized coordinates 
will have been transformed to final form. Thus, much of the 


3. Dynamics of Structures, by W. C. Hurty & M. F. Rubinstein, pp 
ijf. , Prentice Hall, 1964 
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control in obtaining diagonalized dampma and diagonalized stiff- 
ness will have been removed, so it appears that the Generalized 
stiffness matrix may be coupled. Since j^I,] is a square matrix, 
the requirement of equation (23) implies that the order of the 
adioint vectors ft be the same as the order of the primary vectors 

c* 

A. Use Gram-Schmidt method for orthocronalizing the Pri- 
mary Ritz modes and apply the simple constraint of unit diagonal- 
izing with no weighting. 

B. But in the case of Adnoint Ritz modes the orthocronal- 
ized set will be expanded in terms of the complete set of adioint 
bases with the dual constraint of equation (23) havina mass 
weicrhting. 

Mathematically these statements translate into buildincr 
the normalized vectors as follows. 


A. Self -Orthogonal izat ion 


To start with, the simple Gram-Schmidt method sets up a matrix of 
undetermined coefficients to involve an increasing number of 
oases vectors $ in the content of the orthogonal lzed vectors Q. 
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The self normalizing constraint is written in matrix notation as: 

W’kJ ■ [^0 

One can better appreciate the earlier topic of having to decide 
on sequencing during generation of ingredients of the bases in 
light of the character of orthogonalization. If the fundamen- 
tals are sequenced together at first, then the initial orthogo- 
nalized vectors will contain a minimum of higher harmonics in 

their expansion. We ask. Is this good or bad? If inertia 

effects dominate the dynamic behavior of a structure under 

certain loads, probably sequencing harmonics in earlier might 

help. But in this study the option was taken to group the fun- 
damentals ahead of the harmonics instead of layering one fun- 
damental and all of its higher harmonics on top of a second 
fundamental and all of its higher harmonics et cetera. 

Trace the effect of the orthogonality constraint on 
achieving a solution for the undetermined coefficients. Operate 
on the first two equations of the set in equation (24). 

(26) |C 2 } = a n{ <t) i} + {^2} * Im P° se the 
single orthogonality constraint of equation (25). 

(27) 

( 23 ) 


KfK) = s “=> 2 KH-uW * {* 2 }} ■ * 


{♦iKi + M {*2} = a n 




<"2} 
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T , 

{M j ^ 7 1 

Substitute from (26) a^ = (-) ^ - - • The coefficient 

{*1} {*1} 


a^ is now expressed entirely in terms of the set of known norm- 
alized bases vectors C$3. Substitute a^ into the equation above 
and now {£->}■ is known. Turn to the third vector. 




This gives 

3 - . 



cut ing 
vector 


two algebraic equations in the two unknowns a_. ^ and 
these coefficients makes ■{ C j- known. Turn to the 
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‘32) {c i+ i} - + a 2ik} + + + {^i+l} * 


- [kUk, kbkk] 


a 2i 

V 
Si 
( 1 ) 


* C *W a j< 


i+1 


Assemble the i constraint conditions in preparation for sub- 
stituting from equation (32) into them. 

,33) kfkk = kfkk - •••* kfkk - 

which can be combined into 


< 33A ’ [k}?k}Tk}k--k)] kk ■ 


o 

6 

lb) 


Now substitute the expansion in terms of 


,34) [kf-kk .kk [kb kb kbkkJ 


/a, . ' 

a 2i 

Si 

a ii 
l -1 J 


= 0. 


Which can be compressed into 
<35) a jkl =0 - 

But by introducing partitions into t<J>J and £aj it can be written 
more intuitively as 

(36) CC3fk3 1 |l»J 1+1 ][ a )| * 0. = = => 


= 0. 


,T.- 


.T, 


T 

-i *■ r 


= = => HCJjC<D3 i Ca} 1 + CCJ;t0i i+i = 0. = = => CQJ^C<pj.Cai).= < - ) C«t>3 1 + 1 
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from which all of the undetermined coefficients a i ,a 2 ,...,a i can 
be computed and substituted into equation (32) to evaluate 

For solutions of succeeding vectors i+2, i+3, ....n? the 

T 

square coefficient matrix increases incrementally in 

order up to n x n, so that results from the previous calculation 
might be considered for saving. Then the increments can be 
merged into the salvaged core to continue on. Details of the 
strategy will be given in the CODING document. 


B. Dual Orthoqonalization 

The dual orthogonal izat ion of the adjoint bases are organized 
into a full matrix of undetermined coefficients for the expansion 
of normalized vectors £1 as components of the raw adjoint bases i|>. 


(37) 


(CAr,}^ 

™i,n } ; 


*11 *12 *13' 
*21 *22 *23 
*31 *32 *33 


ln 

J 2n 

} 3n 


l*nl *n2 *n3 *nn; 


C *L3 3 

C "tn 3 J 


This can be written in more conventional form as 


38 ) [{^ L i},{ a L 2 },- ‘ -{^Ln}] ' [{^Ll} , * ’ * {^Ln}] 


f *ll *21 , “*nl 
Pl2 f?22 ’ ‘ 'Zn2 
d 13 d 23 ’ ’ * D n3 


*in *2n’ ‘ **nn; 


Note the the b^s in equation <j 8) are the transpose of the 
b^ s in equation ( 37 ). Equation (38) can be condensed to 
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(39) [^Ln] = [^Lp] [^pn] ' in whictx is the transpose of the 

b . . ' s . 
ij 

The mass orthogonality constraint, according to the 

transformation requirements of the dynamic differential equation 

(23), is a relationship between the normalized adjoint and the 

T 

normalized primary bases; i.e. S^. Since the M and C, 

are known, their transposes can be taken immediately, so it 
becomes better strategy to transpose the constraint equation in 
order to obtain & in non-transposed form; i.e. 

,4o> k/KJ Ki] * [V]- 

Now substitute from the expansion equation (39) into the mass 
orthogonality constraint equation (40). Confine to one index at 
a time. Set i = 1. 

(41> [ C Lj] [ M Ll] {^Ll} = [^Lj] [ M LL]Kp]{ P pl} =S jl = { 0 for j * 1 


This will produce a solution to the first column of 0 Next 

set i = 2 and substitute for ^L2* rom equation (39). 

(42) [^Lj] [ M Ll] {^L2} = [ C Lj] KlJK P ]{ 3 P 2} =S j2 = { 0 for j * 2 


This will produce a solution to the second column of 0. A pat- 
tern is now apparent for solving for the complete content of 0 in 


a single operation, by recognizing the 

jVl^ J [^Lp] iS tiie same in ever y equation; only 
of unknown 0's and the columns of the change, 

columns into matrices; i.e. 

««> K 3 ] T KL]Kp][y - [\] • 


Isolate 0 pi - 


coefficient 
the columns 
Combine the 
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1441 M * [\][«£)"lL*Lp] 

Substitute into the defining equation for normalized adioints 

( 29 ) . 

( 45) [ il Ln] = [vJ[ 0 pn] = [\p] [ L \] [ C L] M Ll\ P ] * 

All quantities are now derived for getting a solution of 
the dynamic differential equation by using Ritz modal vectors. 

One thing not taken up however, was the definition of damping so 
as to give as sparse a generalized damping matrix as possible. 

The other topic that still needs addressing is data recovery. 

These topics will be reserved for an extension to the basics as 
developed here. Details of converting this theory to DMAP coding 
has been published in a report entitled "RITZ MODES FOR 
UNSYMMETRIC MATRICES— DMAP CODING OF THE THEORY" by Thomas B. 
Butler. 

The coding was done in 3 steps. (1) Fundamental Primary 
and Adioint modes were obtained from a DMAP ALTER to the Statics 
Rigid Format. The listing of this code is attached as Appendix 
A. >3) Harmonics for the Primary and Adioint sets were coded as 
a pure DMAP approach. The listing of this code is attached as 
Appendix E. (3> Orthogonalization was coded as a pure DMAP 
approach. The listing of this code is attached as Appendix C. 

A simple demonstration problem was used to certify the 
method and the coding. It was run twice with a different 
threshhold value of Kappa to exercise a number of different 
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paths. Details of the demonstration problem are given in Appen- 
dix D. Generalized Mass and Generalized Stiffness that were 
produced in the two different runs are given m Appendix E. One 
run set the linear independence threshhold Kappa to 0.007. Only 
the four fundamentals passed the linear independence test, so the 
generalized mass and stiffness are only of order 4x4. Kappa 
was deliberately set high to a value of 0.95 in the second run so 
as to admit 5 harmonics in addition to the four fundamentals. 
Resulting generalized mass and stiffnesses are of order 9 x 9. 

The generalized mass in both cases was practically unity. 
Off diagonal terms were at least 14 orders of magnitude less than 
those on the diagonal. Marked differences show up in the gener- 
alizeo stiffneses for the two cases. When a mathematically 
logical value of Kappa is used as in the .007 case, the terms on 
the diagonal dominate the off-diagonal terms, implying very weak 
coupling between Ritz modes. This weak coupling could very well 
lustify the use diagonal matrices and so benefit from a decoupled 
solution. When an improbable value of Kappa is used as in the 
0.95 case, off-diagonal terms are large. This observation can 
lead to a manageable criterion for completeness of the modes. 
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CONCLUSIONS 


The method is workable. DMAP codincr has been automated 

4 

to such an extent by using the device of bubble vectors, that 
it is useable for analyses in its present form. This feasibility 
study demonstrates that the Rite Method is so compelling as to 
warrant coding its modules m FORTRAN and organizing the result- 
ing coding into a new Rigid Format. 

Even though this Ritz technique was developed for unsym- 
metric stiffness matrices, it offers advantages to problems with 
symmetric stiffnesses. If used for the symmetric case the solu- 
tion would be simplified to one set of modes, because the adjoint 
would be the same as the primary. Its advantage in eithei type 
of symmetry over a classical eigenvalue modal expansion is that 
information density per Ritz mode is far richer than per eigen- 
value mode; thus far fewer modes would be needed for the same 
accuracy and every mode would actively participate in the re- 
sponse. Considerable economy can be realized in adapting Ritz 
vectors for modal solutions. This new Ritz capability now makes 
NASTRAN even more powerful than before. 


4. Butler, T. G. and t'amidi 
Merging", Proceedings the 


F. R. "Bubble Vector in Automatic 
Fifteenth NASTRAN Colloquium, 19a? 
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APPENDIX A 


DMAP CODING FOR 

PRIMARY AND ADJOINT FUNDAMENTALS 

3$ ADVISABLE TO INCLUDE THIS NASTRAN CARD IN JOBS. 

/*. 

3 NASTRAN MAXFILE3 = 60, FILES = ( INPT, INP1 . INP2 , INP3 , INP4 ) 

3 RITZFUND. DMP $ DMAP ALTER FOR GENERATING FUNDAMENTAL RITZ VECTORS $ 

ALTER 2 $ ALTERS FOR 1988 VERSION OF NASTRAN 

PARAML CONTK/ /^PRESENCE*/ / / /V,N,NOCONTK 3 CONTK IS A DMAP INPUT $ 

$3 $ NOCONTK = -1 IF CONTK IS MISSING 3 

PRTPARM / / 0 / C , N , NOCONTK $DB 

COND ERR0R4, NOCONTK $ ABORTS IF USER OMITS CONTK. 

$3 AFTER GP3 

ALTER 26,26 $ REPLACES PARAM STATEMENT WITH ONE THAT ENABLES MASS GENERATION 

PARAM / / *ADD* I NOMGG/ 1/0 $ ALERTS EMG TO GENERATE MGG 3 

ALTER 39 3 AFTER EMA OUTPUTS KGGX FOR STIFFNESS 

ADD CONTK. KGGX/NSKGG/ 3 THIS IS THE NON-SYMMETRIC STIFFNESS. 

ALTER 45 3 AFTER EMA OUTPUTS MGG 

PURGE MNN , MFF . MAA i NOMGG 3 

ALTER 61,61 3 DON'T PURGE OG. 

PURGE KRR , KLR , QR , DM / REACT / GM / MFCF1 / GO , KOO , LOO . PO , OUUV i OMIT / PS . 

KFS.KSS/ SINGLE 3 

ALTER 62,62 3 ADD MGG TO EOUIV 

EQUIV KGG , KNN / MPCF1 / MGG , MNN / MPCF1 3 

ALTER 65,69 3 REPLACE MCE 2 & SCE1 WITH NON-SYM OPN'S 
VEC USET/GNVEC/*G*/*M*/*N* 3 l'S ON N 

PARTN KGG.GNVEC, /KMM,KNM,KMN,KNNBAR/ -1 $ -1 MEANS THAT GNVEC IS USED 
S3 FOR BOTH ROW AND COL PARTNG , BUT DOES NOT MEAN THAT KGG IS SYMMETRIC 
tlPYAD KNM,GM,KNNBAR/KNNl/0/+l/+l 3 

MPYAD GM.KMN.KNN1/KNN2/+1/+1/+1 3 
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MPY3 GM , KMM , KNN 2 / KNN / 0 $ 

MCE2 USET ,GM, ,MGG, , / KWASTE , MNN , , $ 

LABEL LBL2 $ 

EQUIV KNN, KFF/ SINGLE/ MNN, MFF/ SINGLE $ 

COND LBL3. SINGLE $ 

VEC USET/NFVEC/*N*/*F*/*S* $ l'S ON S 

PARTN KNN,NFVEC , /KFF,KSF,KFS,KSS/-1 $ 

SCE1 USET, ,MNN, , / , , ,MFF, , $ 

ALTER 71,71 $ ADD MFF TO EQUIV 

EQUIV KFF, KAA/ OMIT/ MFF, MAA/ OMIT $ 

ALTER 73,82 $ REPLACE SMP1 .RBMGl ,RBMG2 ,RBMG3 WITH UNSYM OPN'S 
VEC USET/FAVEC/*F*/*A*/*0* $ 

PARTN KFF , FAVEC , /KAABAR ,KOA,KAO ,KOO/ -1 $ 

DECOMP KOO/KOL,KOU/0/0/3,N,KOMIN/S,N,KODET/S,N,KOPWR/S,N, DOSING $ UNSYM 

$$ DECOMP OF KOO 

FBS K0L,K0U,K0A/G0/0/-l $ 

MPYAD KAO, GO, KAABAR/ KAA/ 0/+1/+1 $ 

SMP2 USET, GO, MFF /MAA $ 

LABEL LBL5 $ 

EQUIV KAA, KLL/ REACT/ MAA. MLL/ REACT $ 

COND LBL6, REACT $ 

VEC USET/ALVEC/*A*/*L*/*R* $ 

PARTN KAA , AL VEC , /KLL ,KRL ,KLR ,KRR/ - 1 $ 

RBMGl USET, ,MAA/ . , ,MLL ,MLR ,MRR $ 

LABEL LBL6 $ 

OUTPUT 1 , ,,,,/ /-1/4 $ SET THE DEFAULT LABEL 

OUTPUT 1 MLL, KLL,, ,/ /0 /4 $ FOR ORTHOGONAL I SATI ON AND FINAL CHECK 
DECOMP KLL/LLL .ULL/O/O/S .N.KLMIN/S ,N,KLDET/S ,N,KLPWR/S ,N .KLSING $ 

COND LBL7, REACT $ 

FBS LLL ,ULL .KLR/DM/0/-1 $ 

MPYAD KRL.DM.KRR/X/0/+1/+1 $ 

DECOMP X/KL .XU/O/O/S ,N,XMIN/S ,N,XDET/S .N.XPWR/S ,N,XSING $ 
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DECOMP KRR / KRLL , KRUU / 0 / 0 / 3 , N , KRMIN / S , N , KRDET / S,N,KRPWR/S.N,KRSING $ 

PARAMR / /C,N,DIV/ / / /V,N ,RBEPS/V,N ,XDET/V,N, KRDET $ 

PRTPARM / / O/C ,N ,RBEPS $$ RIGID EODY TRANSFORMATION CHECK 
LABEL LBL7 $ 

3$USER INPUT = CONTK , A matrix of control properties of the same order as KLL. 
$$ = SCVEC, A vector of a leading 1 followed by zeroes of order equal to 

the number of loading subcases. 

= SCADJ , A vector of leading zeroes followed by a trailing 1 of order 
equal to the number of loadina subcases. 

- KAPPA, A real double precision parameter for setting the tolerance of 
the cosine linear independence check. 

= MODSPEC, A parameter for setting the integer maximum of harmonics 
to be generated. 

= LONGONE , A Vector of ones whose length is greater than twice the 
number of subcases times MODSPEC olus one. L= 2(30(1 + MODSPEC) 


$$ 




^>9 
$$ 
$ 3 
$$ 


FBS LLL , ULL , MLL / SOLP $ INERTIA COEFFICIENT FOR HARMONIC RITZ P VECTORS 

$$ PER EQUATION (15) 

$ INTRODUCE THE SOL'N OF THE TRANSPOSE OF KLL TO DEVELOP INVERSE VECTORS 
TRNSP KLL/KLLT/ $ 

DECOMP KLLT / TLLL , TULL / 0 / 0 $ 

FBS TLLL, TULL, MLL /SOLA $ INERTIA COEFFICIENT FOR HARMONIC RITZ A VECTORS 

$$ PER EQUATION (15) 

ALTER 88,92 $ REPLACE SSG3 WITH NON SYM OPN'S 
FBS LLL , ULL , PL / ULV / $ 

FBS TLLL, TULL, PL /TULV/ $FAY THE PRICE OF DOING FBS ON INVERSE TOO 

COND NOUOOV,OMIT $ 

FBS KOL ,KOU,PO/UOOV/ $ 

LABEL NOUOOV $ 
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$ 

$$ CALCULATION OF RULV IN CONNECTION WITH IRES WILL BE IGNORED BECUZ IT WOULD 
$$ BE ADVISABLE TO DO THIS CHECKING WITH SYMMETRIC MATRICES. 

$ ======> USER SHOULD ESCHEW USE OF IRES <======= 

$ THE EPSILON SUB E CHECK WILL BE DONE IN DBL PREC BECUZ PAR AMD IS NOW AVAILABLE 
$ 

MPYAD KLL , ULV , PL / DELPL /0/-1/+1 $ 

MPYAD ULV, DELPL, /DELWORK/+1/+1 $ 

MPYAD PL, ULV, /ALLWORK/ +1/+1 $ 

SCALAR DELWORK/ / 1/ 1/ /V,N,EPSNUM $ 

SCALAR ALL WORK/ / 1/ 1 / /V.N.EPSDEN $ 

PAR AMD / / *DI V* / V , N , EPSUBE / V , N , EPSNUM / V , N , EPSDEN $ 

PRTPARM / /0/C, N, EPSUBE $$ RIGID BODY TRANSFORMATION CHECK 

PARAM / /*ADD*/V,N,ADJCYC/+l/0 $ VALUE OF PARAM REMAINS POSITIVE DURING 

3$ PROCESSING OF PRIMARIES. 

COPY LONGONE/CLONONE/ 0 $ 

ADD LONGONE , CLONONE / LONGNULL / ( - 1 . 0 , 0 . 0 ) $ 

PARAML SCVEC / /^TRAILER*/ 2/V.N.VECRO $ ROW SIZE IS READ FROM 3CVEC 

PARAML ULV/ / ^TRAILER*/ 1 / V ,N ,ZCOL $ COL SIZE IS READ FROM ULV 

PARAM / / *EO * / V , N , LODNO / V , N , VECRO / V , N , ZCOL $ LODNO IS NEGATIVE IF VECRO=ZCOL 

COND ADJLUP, LODNO $ CONTINUE IF SCVEC AND ULV AGREE 

JUMP ERR0R4 $ ABORT IF SCVEC AND ULV DON'T AGREE 

FILE FALCCLI= SAVE/FALRRLI= SAVE $ 

LABEL ADJLUP $ TOP OF LOOP PRIOR TO FORMATION OF EITHER PRIM OR ADJ FUND 

COND AD JTRN , ADJ CYC $ ADJCYC INITIALLY IS POS TO GIVE PRIORITY TO PRIMARY 

COPY ULV/FLV/ 0 $ FLV WILL REMAIN AN INTERNAL DATA BLOCK TO THIS FUND LOOP 

EQUIV ULV, CLONFLV $ CLONFLV IS INTERNAL. EQUIV WILL BE BROKEN AT TOP'O LOOP 

JUMP PRIMSEG $ GO AROUND THE AJOINT PREP 
LABEL ADJTRN $ 

COPY TULV/FLV/ 0 $ FLV IS INTERNAL 

EQUIV TULV, CLONFLV $ CLONFLV IS INTERNAL 

LABEL PRIMSEG $ 
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MPYAD CLONFLV,FLV r /FLSQ/+1 $ 

DIAGONAL FLSQ/SCALF/ *DI AGONAL* /C ,N, -0.5 $ 

MPYAD FLV.SCALF, /FMOD/ 0 $ CANDIDATE MODE NORMALIZED TO EUCLID LGTH. 

COFY FMOD/FCLON/ 0 $ 

MPYAD FCLON.FMOD, /FD0T/+1 $ MATRIX ORDERS = ZCOL. 

MATPRN FDOT .,,,11 $DB 

COPY SCVEC / SCVEC I / 0 $ ROW BUBBLE STARTING FROM HEAD 

COPY SCVEC / LI VECI / 0 $ COL BUBBLE STARTING FROM HEAD 

COPY SCVEC /MODPARTN/ 0 $ DUMMY TO BE USED FOR SWITCHING WITHIN LOOPS 

ADD LIVECI , SCVEC /FALCCL I / /(-l. 0,0.0) $ NULL BUT SAME LENGTH AS SCVEC 

COPY FALCCL I / FALRRL I / -1 $ NULL SAME LGTH AS SCVEC 

PARAM / /*MPY*/V,N,ROCNT/ 1 / 0 $ RESET ROW COUNT TO 0 BEFORE LI CHECK 

LABEL LIRLUP $ TOP OF ROW PORTION OF LINEAR INDEPENDENCE LOOP 

PARAM . / /*MPY*/V,N,NORFAL/+l/-l $ SET DEFAULT TO NEG TO JUMP OVER FAIL BOOK 

PARTN FDOT, ,SCVECI/ .ROCAI, , /+7/+2 $ 

PARAM / / *ADD*/V,N,ROCNT/V,N,ROCNT/ 1 $ ROW COUNT MONITOR INCREMENTED BY ONE 

PARAM / / *MPY*/V ,N,CLCNT/ i /V,N,ROCNT $SET COL COUNT=ROW COUNT PRIOR LI CHK 

LABEL LICLUP $ TOP OF COLUMN PORTION OF LINEAR INDEPENDENCE LOOP 

PARAM / / *ADD*/V,N,CLCNT/V,N,CLCNT/ 1 $ COL COUNT MONITOR INCREMENTED BY ONE 

SCALAR ROCAI/ / 1/V,N,CLCNT/ /V,N,RCF $ COSINE TERM TO BE TESTED 

PARAMD / /*ABS*/V,N.COSRCF/V,N,RCF $GETS ABSOLUTE VALU OF COS (ROW, COL) TERM 

PARAMD //*LE*/ /V,Y,KAPPA/V.N,COSRCF////V,N,LICHK $LICHK =-l IF KAPPA < COSRCi 

PARTN LIVECI. .SCADJ/CDUM, , ,/+7/+l $ HAVE BUBL VEC TO TRACK. TRIM TRAIL ZER( 

MERGE CDUM, , , , , SCVEC / L I VEC J / + 7 / + 1 $BUBBLE INCREMENTED AWAY FROM HEAD 

COND FALBOOK , L ICHK $ CATALOG FAILURE POSITION 

JUMP MORCLI $ SKIP AROUND CATALOGING IF TEST WAS SUCCESSFUL 

LABEL FALBOOK $ 

PARAM / /*MPY*/V,N,NORFAL/V,N,LICHK/ -1 $ SETS SIGNAL ONLY WHEN A COL FAILS. 
$$ HAS OPPOSITE SIGN TO LICHK. POSSIBLE REPEATS ARE O.K. 

PRTPARM / /0/C,N,R0CNT $ ROW # OF CANDIDATE WHICH FAILED LI TEST 
PRTPARM / / 0 / C , N , CLCNT $ COL # OF CANDIDATE WHICH FAILED LI TEST 
ADD FALCCLI , LIVECJ/ FALCCL J / $ ACCUM OF COL POS'NS OF FAILURES 
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SWITCH 

LABEL 

SWITCH 

PARAM 

COND 

REPT 

■P 

LABEL 

A, 

9 

COND 
ADD 
SWITCH 
LAEEL 
PARAM 
PARAM 
COND 
JUMP 
LABEL 
PAR AML 
PAR AML 
PARAM 
$$ 

PRTPARM 

COND 

SWITCH 

JUMP 

LABEL 

SWITCH 

LABEL 

PARTN 

PARTN 

$$ 


FALCCLI , FALCCLJ / / V,N,LICHK $ 

MORCLI $ CONTINUE LI TESTING IN THIS ROW EVEN AFTER A COL FAILS 
LIVECI ,LIVECJ/ / -1 $ 

/ /*EQ*/V.N,LICDUN/V.N,CLCNT/V,N,ZCOL $ LICDUN = -I IF CLCNT = ZCOL 
NUROW.LICDUN $ JUMP OUTSIDE OF COL LOOP IF (3 LAST COL 
LICLUP ,999 $$$$$$$$$$$ END OF COLUMN LOOP 

NUROW $ 

GUDRO.NORFAL $ JUMP IF NO COLS IN CURRENT ROW HAD A LI FAILURE 
FALRRLI, 3CVEC I / FALRRL J / $ ACCUMULATED ROW POSN'S OF FAILED ROWS 
FALRRLI , FALRRLJ / / -1 $ SWITCH ONLY IF THIS ROW FAILED, ELSE STAYS 
GUDRO $ 

/ /*SUB*/V.N,ROTEST/V,N.ZCOL/i $ DECREMENT ZCOL BY ONE FOR ROTEST 
/ / *EQ* / V , N , L IRDUN / V , N , ROCNT / V , N , ROTEST $ LIRDUN = -l.IF R0CNT=R0TE3T 
KLENUP , L IRDUN $ 

MOROW $ GO AROUND CLEAN UP IF MORE ROWS REMAIN TO BE TESTED 
KLENUP $ 

FALCCLI/ / ^TRAILER* /C,N,6/V,N, FALCDENS $ DENSITY OF THE COL FAIL VEC 
FALRRLI/ / *TRAI LER* / C , N , 6 / V , N , FALRDEN S $ DENSITY OF THE ROW FAIL VEC 
/ /*LE*/V.N,DENSLK/V,N,FALCDENS/V,N,FALRDENS $ IF COL DENS (/= ROW DEN. 

DENSITY SELECTION PARAMETER IS NEGATIVE 
/ / C , N , 0 / C , N , DEN3LK $DB 

OTHER, DENSLK $ SWITCH MODPARTN TO THE VECTOR W LWR DENSITY 
FALRRLI ,MODPARTN/ /-I $ 

MODSET $ 

OTHER $ 

FALCCLI, MODPARTN/ /-I $ 

MODSET $ 

FMOD, MODPARTN, /PHI, , , /+7/+2 $SURVIVORS OF LI TEST 
SCALF, , MODPARTN / PFVEC , , , /+7/+1 $ CLUSTER VECTOR FOR MERGING FUND- 
AMENTAL MODES WITH HARMONICS 
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MERGE PFVEC, , , , ,LONGNULL/HEADPF/+7/+l $PUT PFVEC AT THE HEAD OF A LONG VEC 
PARTN LONGONE. , HEADPF / SHORTONE , . , /+7/+1 $ PARTITION LONGONE DOWN 

PERMANENT 

MERGE SHORTONE, , , , .LONGNULL/HEDSHORT/+7/+1 $ APPEND PERMANENT-SIZE NULL TC 
TAIL 

ADD HEDSHORT, LONGONE /NEGTAL/ (1.0,0.0)/(-l. 0,0.0) $PERMANENT-SIZE NEG 0 TAII 

$ 

$ PROVIDE FOR THE POSSIBILITY OF THE CULLING VECTOR CONTAINING l'S IN THE END 
$ POSITIONS, WHICH WOULD DESTROY THE FUNCTION OF THE SHIFTING VECTORS. CONVERT 
$ SCALF TO ALL ONES. 

A* 

■9 

MERGE SCALF SCADJ/COLSCAL/+7/-)-2/+2 $ SETS TRAILER TO RECTANGULAR 

TRNSP COLSCAL/SCALFRO $ CONVERT COL TO ROW 

MPYAD COLSCAL , SCALPRO , /SQUID/0 $SOUARE MTX OF READ D.P. IN PREP FOR DIAG 

DIAGONAL SOUID/ FULU / /'COLUMN* / 0 . 0 5FULU IS A CLUSTER OF ALL I '3. LGTH=FMOD 
PARTN FULU, ,MODPARTN/FUNPART. TOSS, ./+7/+1/+2/+2 $LGTH 1ST=PHI ,T0SS=C0MP WR" 

FMOD 

MERGE FUNPART .FULU/HEDCLUS/+7/+1 $FORM CLUSTER OF FUNPART 0 VEC HEAD 

MERGE TOSS. , , , .FULU/HEDTOSS/+7/+1/+2 $FORM CLUSTER OF TOSS Q VEC HEAD 
ADD FULU, HEDT0SS/TALCLUS/(i.0,0.0)/(-i. 0,0.0) $ CAP ZEROES 0 HEAD Of 

FUNPART 
$ 

COND AD JWR AP , AD J CYC $ 

$ 

COPY PHI/PHIPI/ 0 $ 

ADD NEGTAL. / LONGPRMI / ( - 1 . 0 , 0 . 0 ) $ TAIL CLUSTER = PERMANENT PRIMARY MODES 

PARTN SCVEC , .HEDCLUS/ .HMHED, . /+7/ -hi $ HEAD SHIFTER TRIMMED TO LGTH = PHIPI 

PARTN SCADJ , .TALCLUS/ ,HMTAL, , /+7/+1 $ TAIL SHIFTER TRIMMED TO LGTH = PHIPI 

JUMP ADJKUNT $ 

LABEL ADJWRAP $ 

COPY PHI/PHIAI/ 0 $ 

ADD NEGTAL, /LONGARMI/ ( -1. 0,0. 0) $ TAIL CLUSTER = PERMANENT ADJOINT MODES 
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PARTN SCVEC , .HEDCLUS/ .ADJHED, , /+7/+1 $ SHIPPER TRIMMED TO LGTH = PHIAI 

PARTN SCADJ , , HEDCLUS / , ADJTAL , , / +7 / +• 1 $ SHIFTER TRIMMED TO LGTH = PHIAI 

$ 

JUMP HARMONY $ 

A 

LABEL MOROW $ 

PARAM / /*MPY*/V,N,CLCNT/l/0 $ RESET COLUMN COUNT TO ZERO 

PARTN SCVEC I, .SCADJ/RDUM, , , /+7/+1 $ TRIM TRAILING ZERO 

MERGE RDUM, , , , .SCVEC/SCVECJ/+7/ +1 $ BUBBLE INCREMENTED AWAY FROM HEAD 

SWITCH SCVEC I , SCVEC J / / -1 $ 

COPY 5CVECI /LIVECI / 0 $ COL BUBBLE INDEX ALIGNED WITH ROW TRACKER 

$ 

REPT LIRLUP, 999 $ 

$ 

i LABEL ADJKUNT $ 

! PARAM / /*MPY*/V,N,ADJCYC/1/-1 $ 

$ 

REPT AD JLU P.999 $ 

A 

LABEL HARMONY $ 

3 

OUTPUT 1 1 /-I/O $ CALLS THE DEFAULT LABEL. NEEDED FOR REWINDS LATER. 

OUTPUT! PHIPI .PHIAI, SOLP , SOLA, / /O/O $MANY CALLS TO BE MADE IN HARMONIC PHASE 

OUTPUTI I 1-1 / 1 $ SETS THE DEFAULT LABEL 

OUTPUT1 HMHED, HMTAL, AD JHED, ADJTAL, / /O/l $ 

OUTFUT1 , ,,,,/ / -1/2 $ SETS THE DEFAULT LABEL 
OUTPUTI LONGNULL , LONGPRMI , LONGARMI , LONGONE , / /0/2 $ 

S 

ALTER 154,154 $ REMOVE OPTIMIZATION LOOP TO PREVENT 0 ' FLOW OF CEITBL 
$ 

ENDALTER $ END OF ALTER FOR RITZ FUNDAMENTAL MODES 
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APPENDIX B 


DMAP CODING FOR 
PRIMARY AND ADJOINT HARMONICS 

AFP DMAP $ FOR EXECUTION AFTER RITZFUND. INPUTS FROM INPT, INP1 , INP2 
BEGIN $ PROGRAMMED FOR 1988 VERSION OF NASTRAN. OUTPUT TO INP3 
FILE LONGPRMI = SAVE / LONG ARM I = SAVE / PGEN I = SAVE / AGEN I = SAVE / FULU = SAVE $ 

FILE PHHED = SAVE / PHTAL = SAVE / AHHED= S AVE / AHHED = SAVE / AHTAL = SAVE $ 

PARAM / /*MPY*/V,N,FRIMCYC/+l/-i $ CONTROL PARAM FOR PRIMARY 1ST HARMONIC 
PARAM / /*ADD*/V.N,ADJCYC/+1/ 0 $ CONTROL PARAM FOR ADJUNCT 1ST HARMONIC 

PARAM / / * ADD* / V , N , NUPGEN / +1 / 0 $ CONTROL PARAM FOR PRIMARY HIGHER HARMONICS 

PARAM / /*ADD*/V,N,NUAGEN/+1/ 0 $ CONTROL PARAM FOR ADJUNCT HIGHER HARMONICS 

PARAM / / *MP Y */V,N, HARMNO / 1/ 0 $ SET THE HARMONIC COUNTER TO ZERO. 

$ 

LABEL HMNICGEN $ TOP OF LOOP FOR HARMONIC GENERATION %%%%%%%%%%%%% 

$ 

PARAM i i *ADD* / V , N , HARMNO / V , N , HARMNO / 1 $ INCREMENT THE HARMONIC COUNT BY ONE 
COND PRIMPREP .PRIMCYC $ 

COND ADJPREP , ADJCYC $ 

COND PHMNPREP, NUPGEN $ 

COND AHMNFREP , NUAGEN $ 

LABEL PRIMPREP $ 

INPUTT1 /. , . , /-I/O $ 

INPUTT1 /HLV /O/O $ READ PHIPI INTO HLV 

COPY HLV/ PHIPI/ 0 $ THIS IS THE ROOT FOR THE 1ST MERGE OF HARMONICS TO PRIM! 

EQUIV HLV, TESTER /PRIMCYC $ HLV AND TESTER ARE INTERNAL NAMES OF GENERATOR. 

$$ LATER ON, EQUIV WILL BE BROKEN AT TOP'O LOOP. 

INPUTT1 /KMMTX, , , , /I/O $ SKIP PASSED 2ND DB AND READ SOLP INTO KMMTX 
INPUTT1 /, , , , /-l/l $ 

INPUTT1 /HEDVECI ,TALVECI , , , /0/1 $READ HMHED INTO HEDVECI & HMTAL INTO TALVECI 
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INPUTTI I ... . /-l/l $ 

INPUTT 1 / PGENI , PHTAL , , , / 0 / 1 $READ HMHED INTO PGENI ;READ HMTAL INTO PHTAL. 

COPY PGENI /PHHED/ 0 $ DUMMY STATEMENT TO FOOL THE COMPILER 

ADD HEDVEC I , TAL VEC I / FULU $ FLUFF TO HELP FIAT LOCATE THE REAL FULU 

JUMP HMYBUS $ GO AROUND THE A JOINT PREP 

LABEL ADJPREP $ 

INPUTTI -I/O $ REWIND FROM PREVIOUS PASS THRU LOOP AND POSITION @ 1ST DB 

INPUTTI /HLV, , , , /l/O $SKIP PASSED 1ST DB AND READ PHIAI INTO HLV 
COPY HLV/ PHIAI /0 $THIS IS THE ROOT FOR THE 1ST MERGE OF HARMONICS TO ADJOIN: 
EQUIV HLV, TESTER/ AD JCYC $ HLV & TESTER ARE INTERNAL NAMES OF GENERATOR 
INPUTTI / KMMTX , , , , /l/O $ SKIP PASSED 3RD DB AND READ SOLA INTO KMMTX 
INPUTTI /, , , , /-l/l $ 

INPUTTI /HEDVECI ,TALVECI , , , / 2 / 1 3 READ ADJHED INTO HEDVECI & ADJTAL INTO TAL VEC. 
INPUTTI , /-l/l $ 

INPUTTI / AGENI , AHTAL , , , / 2 / 1 SSKIP 2DB & READ ADJHED > AGENI ;READ ADJTAL > PHTAL . 

COPY AGENI /AHHED/ V.N.ADJCYC $ 

JUMP HMYBUS $ 

LABEL PHMNPREP $ 

COPY PGENI /HLV/ 0 $ 

EQUIV PGENI, TESTER /NUPCEN $ TESTER IS INTERNAL NAME OF GENERATOR 
INPUTTI / .... /-I/O $ REWIND FROM PREVIOUS PASS THRU LOOP AND POSITION @ 1ST DI 
INPUTTI /KMMTX,,,, /2/0 $ SKIP PASSED 1ST 2 DB'S AND READ SOLP INTO KMMTX 
EQUIV PHHED, HEDVEC I /NUPGEN/ PHTAL, TAL VEC I /NUPGEN $ 

JUMP HMYBUS $ 

LABEL AHMNPREP $ 

COPY AGENI /HLV/ 0 $ 

EQUIV AGENI, TESTER /NUAGEN $ TESTER IS INTERNAL NAME OF GENERATOR 
INPUTTI / .... /-I/O $ REWIND FROM PREVIOUS PASS THRU LOOP AND POSITION 0 1ST DI 
INPUTTI /KMMTX,,,, /3/0 $ SKIP PASSED 1ST 3 DB'S AND READ SOLA INTO KMMTX 
EQUIV AHHED , HEDVEC I / NUAGEN / AHTAL , TAL VEC I / NUAGEN $ 

LAEEL HMYBUS $ 
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PARAM / / *MPY*/ V,N , INITIAL/ +1 / - 1 $ NEWLY GENERATED CANDIDATE GOING TO 
$$ INITIAL TEST. SET TO -1. 

MPYAD KMMTX,HLV, /UHC/0/1/0/2 $ CANDIDATES OF THE INERTIAL RESPONSES 
PURGE KMMTX $ 

COPY UHC / CLONUHC / 0 $ 

MPYAD CLONUHC, UHC, /UHCSQ/ +1/+1/0/2 $ 

DIAGONAL UHCSQ /SCALHMDIAGONAL*/- 0.5 $ VECTOR OF EUCLIDEAN LENGTHS. 

MPYAD UHC , SCAE.H , /PHIHC/0/1/0/2 $ MATRIX OF CANDIDATE FIRST HARMONIC MODES. 

COPY PHIHC/CANDIDAT/ 0 $ USE GENERALIZED LOOP NAMES 

$$ 

LABEL LIPREP $ TOP OF LOOP FOR LINEAR INDEPENDENCE CHECKING %%%%%%%%%%%%% 
PARAM / /*MPY*/V,N,RFALNO/l/-l $ SET DEFAULT TO NEGATIVE. 

A A 

MPYAD CAND I DAT, TESTER, /CVST/ + 1/+JL/0/2 $ HARMONICS AGAINST THE GENERATORS. 

MATPRN CVST,,,,/ / $ 

COND RECTO .INITIAL $ SUBSTITUTE A REDUCED FULU IF REMNANT < GENERATOR 

DIAGONAL CVST/ FULU /^COLUMN*/ 0 . 0 $ALL ONE VEC FOR HARM VS HARM (CVST IS SQUARE) 
LABEL RECTO $ 

EQUIV HLV, TESTER /MODSPEC $ BREAK EQUIV W TESTER 

EQUIV PGENI .TESTER/ MODSPEC /AGEN I .TESTER /MODSPEC $ BREAK EQUIV W TESTER 

COPY HEDVEC I / HMROWI / 0 $ ROW BUBBLE STARTING FROM HEAD 

ADD HMROWI , HEDVEC I / FALHR 1/ /(-l. 0,0.0) $ NULL . SAME LENGTH AS HEDVEC I 

COPY FALHR I / FALHCI / 0 $ NULL. SAME LENGTH AS HEDVECI 

PARAM / /*MPY*/V,N,CLKNT/ 1/ 0 $ RESET COL COUNT TO 0 BEFORE LI CHECK 

PARAM / /*MPY*/V,N,ROKNT/ 1/ 0 $ RESET ROW COUNT TO 0 BEFORE LI CHECK 

PAR AML CVST//*TRAILER*/2/V,N,HROW $ ROWS IN CVST 

COND RECTT, INITIAL $ JUMP IF ON RECTANGULAR CYCLE 

PARAM / /*SUB*/V,N,HROW/V,N,HROW/l $ REDUCE ROW TEST VALUE 1 FOR TRIANGLE 
LABEL RECTT $ 

PAR AML CVST//*TRAILER*/l/V,N,HCOL $ COLS IN CVST 
$ 

LAEEL HLIRLUP $ TOP OF ROW PORTION OF LINEAR INDEPENDENCE CHECK **%%%%%%*%% 
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COPY HEDVEC I / HMCOL I / 0 $ COL BUBBLE STARTING FROM HEAD 

PARTN CVST, ,HMROWI / ,ROCH, , Z+7/+2 $ 

PARAM / /*ADD*/V,N,ROKNT/V,N,ROKNT/ 1 $ ROW COUNT MONITOR INCREMENTED BY ONE 
COND RECTI, INITIAL $ 

PARAM / /*MPY*/V,N,CLKNT/l/V,N,ROKNT $ SET COL COUNT=ROW COUNT IF TRIANGLE 
LABEL RECTI $ 


LABEL 


HLICLUP $ TOP OF COLUMN PORTION OF LI LOOP %%%%%%%%%%%%%%%%%%%%%%%%%% 


9 

PARAM / / * ADD* / V , N , CLKNT / V , N , CLKNT / 1 $ INCREMENT THE COLUMN COUNT 
SCALAR ROCH/ / 1/V,N, CLKNT/ /V,N,RCH $ RCH IS DBL PREC FETCH OF COSINE 
$$ TERM IN ROW POSITION INDEXED BY THE CONSTANT PARAMETER '1' AND IN THE COL 
$$ POSITION INDEXED BY THE VARIABLE PARAMETER CLKNT. 


PAR AMD 
PARAMD 
COND 
JUMP 
LABEL 
COND 
PARAM 
$$ 

LABEL 

PRTPARM 

PRTPARM 

COND 

ADD 

SWITCH 

$ 


/ /*ABS*/V,N,COSRCH/V,N,RCH $ GETS ABS . VAL. OF COS(ROW,COL) TERM 
/ /*LE*/ /V,Y, KAPPA/ V,N,COSRCH/ / / /V,N, LIHZK $ 

CATALOG, LIHZK $ GO TO CATALOGING IF LIHZK IS NEGATIVE 
MORHCOL $ JUMP TO MORE COL PROCESSING IF TEST PASSED 
CATALOG $ 

RECT2 , INITIAL $ 

/ / *MPY * / V , N , RFALNO / V , N , LIHZK / -1 $ SETS SIGNAL ONLY WHEN A COL FAILS. 

RFALNO TAKES ON OPPOSITE SIGN TO LIHZK. REPEATS ARE OK. 

RECT2 $ 

/ / O/C ,N ,ROKNT $ ROW # OF CANDIDATE TERM WHICH FAILED LI TEST 
/ /0/C,N, CLKNT $ COL # OF CANDIDATE TERM WHICH FAILED LI TEST 
RECT3, INITIAL $ 

FALHCI ,HMCOLI/FALHCJ / $ ACCUMULATION OF COL POS'NS OF FAILURES 
FALHCI ,FALHCJ / /V,N, LIHZK $ 


349 
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LABEL 

$ 

PAR AM 
COND 
PARTN 
MERGE 
SWITCH 
REPT 
$ 

LABEL 

$ 

COND 

LABEL 

ADD 

SWITCH 

COND 

PARAM 

LABEL 

$ 

LABEL 

$ 

PARAM 

COND 

PARAM 

PARTN 

MERGE 

SWITCH 

$ 

REPT 

$ 


MORHCOL $ CONTINUE LI TESTING IN THIS ROW EVEN AFTER A COL FAILS 

/ /AEQA/v,N,LICDON/V,N,CLKNT/V,N,HCOL $ LICDON = -1 IF CLKNT = HCOL 
GNROW,LICDON $ JUMP OUTSIDE OF COL LOOP IF 0 LAST COL. 

HMCOLI , , TALVECI /DUMMY, , ,/+7/+l $ TRIM TRAILING ZERO 

DUMMY, , , , ,HEDVECI/HMCOLJ/+7/+l $ BUBBLE INCREMENTED AWAY FROM HEAD 
HMCOLI ,HMCOLJ / / -1 $ 

HLICLUP , 999 $ $$$$$$$$$$$$$$$ END OF COLUMN LOOP 
GNROW $ CONSIDER TESTING ANOTHER ROW 

GODRO ,RFALNO $ JUMP IF NO COLS IN CURRENT ROW HAD A LI FAILURE 
RECT3 $ CATALOG ROW FAILURE 

FALHRI, HMROWI /FALHRJ / $ ACCUMULATED ROW POS ' NS OF FAILED ROWS 
FALHRI , FALHRJ / / -1 $ SWITCH ONLY IF THIS ROW FAILED, ELSE STAYS 
RECT4, INITIAL $ BYPASS IF ON RECTANGULAR ROUTE 
/ /AMPYA/V,N,RFALNO/V,N,RFALNO/-L $ RESET TO NEGATIVE. 

RECT4 $ 

GODRO $ 

/ / AEQ a / V , N , L IRDON / V , N , ROKNT / V , N , HRO W $ LIRDON = -1 IF ROKNT = HROW 

CLENUP , LIRDON $ JUMP OUT OF LI CHECKING IF MTX IS COMPLETELY EXAMINED 

/ /AMPYA/V,N,CLKNT/l/0 $ RESET COLUMN COUNT TO ZERO 

HMROWI, , TALVECI /DUMMY, , , /+7/+1 $ TRIM TRAILING ZERO 

DUMMY, ,,, ,HEDVECI/HMR0WJ/+7/+l $ BUBBLE INCREMENTED AWAY FROM HEAD 

HMROWI, HMROWJ/ /-I $ 

HLIRLUP ,999 $ END OF ROW PORTION OF LINEAR INDEPENDENCE LOOP %%%%%%%%%; 


350 


RITZ MODES FOR UNSYMMETRIC MATRICES 
DMAP CODING OF THE THEORY 


LABEL CLENUP $ 

$ 

COND RECT5 , INITIAL $ 

PARAML FALHCI/ /*TRAILER*/C,N,6/V,N,DENSFALC $ DENSITY OF THE COL FAIL VEC 

PARAML FALHRI / /*TRAILER*/C,N,6/V,N,DENSFALR $ DENSITY OF THE ROW FAIL VEC 

PARAM / / *LE* / V ,N , SLKDENS / V , N , DENSFALC / V , N , DENSFALR $ IF COL DENS 

$$ </= ROW DENS, THE DENSITY SELECTION PARAMETER IS NEGATIVE. 

COND UTHER, SLKDENS $SET MODEPART TO THE VECTOR WITH LWR DENSITY 
LABEL RECT5 $ 

COPY F ALHR I / MODEF ART / 0 $ 

JUMP MODESET $ 

LABEL UTHER $ 

COPY FALHCI /MODEPART/ 0 $ 

LABEL MODESET $ 

PARAML MODEPART/ /*TRAILER*/C,N,6/V,N,M0DENSY $ DENSITY OF MODEPART 
PRTPARM / / 0 / C , N , MODENS Y $ 

PARAM / /*GE*/V,N,FILLED/V,N,MODENSY/10000 $ FILLED*- 1 IF MODEPART IS FULL 
COND FOLD, FILLED $ 

JUMP FLEDGE $ 

LABEL FOLD $ SAVE AND GO 

PARAM / / *GT* / V , N , 3UMHUM/ V , N , HARMNO / 2 $ SUMHUM=-1 IF 1ST HARMS OF P £ * 

PASSED 

PRTPARM / /0/C,N,SUMHUM $ 

COND ORTHOG , SUMHUM $ IF 1ST HARMS FAIL RESTORE ORIGINAL NAMES TO OUTPUT 

INPUTT1 /, , , , /-l/O $ 

INPUTT1 /PHIPI ,PHIAI , , , / 0 / 0 $ FUNDS W/O HARMONICS 
INPUTT1 /, , , , / -1/1 $ 

INPUTT1 /HEADPHI ,TAILPHI , , , / 0 / 1 $ TRACKERS W/O HMNIC . HMHED=HEADPHI . HMTAL=TAILPHj 
JUMP N0BI2NEZ $ COPY OUT AS THEY CAME IN 
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LABEL FLEDGE $ 

PARTN C AND I DAT, MODEP ART, /REMNANT, , , /+7/+2 $SURVIVORS OF LI TEST 

$$ FROVIDE FOR THE POSSIBILITY OF THE CULLING VECTOR CONTAINING 1 ' S IN THE 
$$ ENDS, WHICH WOULD DESTROY THE SHIFTING VECTORS. CONVERT SCALH TO ALL ONES. 
COND RECT6, INITIAL $ 

JUMP TRIA6 $ 

LABEL RECT6 $ 

MERGE SCALH, ,,, ,TALVECI/COLSCAL/+7/+2/+2 $ SETS TRAILER TO RECTANGULAR 
TRNSP COLSCAL/SCALPRO $ CONVERT COL TO ROW 

MPYAD COLSCAL , SCALPRO , / SQUID/ 0 $ SQAURE MTX OF REAL S.P. IN PREP FOR DIAG 
DIAGONAL SQUID / FULU / ^COLUMN* / 0 . 0 $FULU IS A CLUSTER OF ALL ONES LGTH=CANDIDAT 
LABEL TRIA6 $ 

PARTN FULU, ,MODEPART/HMYPART,TOSS, ,/7/l/2/2 $1'3.LGTH 1ST =REMNANT.TOSS=COMFI 
MERGE HMYPART, , , , ,FULU/HEDCLUS/ +7 /+1 / 2 $FORM CLUSTER OF HMYPART 0 VEC HEAD 
PARTN HEDVECI , ,HEDCLUS/ ,HEDVECJ ,,/7/l/2/2 $NU SHIFTER HAS LGTH = REMNANT 
MERGE TOSS, ,,, ,FULU/HEDTOSS/ 7/1/2 $FORM CLUSTER OF TOSS 0 VEC HEAD 

ADD FULU, HEDTOSS /TALCLUS/ ( 1. 0,0. 0) /( -1.0,0. 0) $CAP ZEROES 0 HEAD OF HMYPAR: 

PARTN TALVECI , , TALCLUS / ,TAL'7ECJ , , / 7 / 1 / 2 / 2 SNU SHIFTER HAS LGTH = REMNANT 

SWITCH HEDVECI ,HEDVECJ / / -1 $ 

SWITCH TALVECI .TALVECJ / / -1 $ 

PRTPARM / /0/C,N, INITIAL $ 

COND TEST2, INITIAL $ If INITIAL IS NEGATIVE GO TO 2nd LI TEST 

JUMP DEPOT $ READY FOR MERGING AND GENERATING (l.ADJ 2 . NU PHMNY 3.NU AHMNY ) 

LABEL TEST2 $ 

$Test whether one or less columns of REMNANT are left. SET PARAMETER IF SO. 

PAR AML REMNANT/ /*TRAILER*/l/V,N,REMCOL $ 

PARAM / / ALE* / V , N , SALVAGE / V , N , REMCOL / 1 $ 

COND DEPOT, SALVAGE $ 

SWITCH REMNANT, CANDIDAT/ /V,N, INITIAL $ 

COPY CANDIDAT/TESTER/V,N, INITIAL $ 

PARAM / /*MPY*/V,N,INITIAL/V,N, INITIAL/ -1 $ 

REPT LIPREP,! $ MAKE A 2ND PASS FOR THE LI TESTS ON HARMONICS ALONE. 
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$ 

LABEL 

COND 

JUMP 

LABEL 

INPUTT1 

INPUTT1 

EQUIV 

JUMP 

LABEL 

COND 

JUMP 

LABEL 

INPUTT1 

INPUTT1 

EQUIV 

JUMP 

LABEL 

COND 

JUMP 

LABEL 

EQUIV 

JUMP 

LABEL 

COND 

JUMP 

LABEL 

EQUIV 

LABEL 

INPUTT 1 

INPUTT1 

INPUTT1 


DEPOT $ STAGING POINT. WRAPUP. MERGE. HARMONIC GENERATION. 
QUIV1 ,PRIMCYC $ 

CHKADJ $ 

QUIV1 $ 

I , , , , 1-112 $ 

/LONGPRMI , , , , / 1/2 $SKIP 1 DB & READ IN LONGPRMI 
LONGPRMI , LONGRMI / PR IMC YC / PH IP I , PH I I / PR IMCYC $ 

MERGBUS $ 

CHKADJ $ 

QUIV2 , ADJCYC $ 

CHKPGEN $ 

QUIV2 $ 

/, , , , / - 1 •' 2 $ 

/ LONGARMI , , , , 1212 $SKIP 2 DB & READ IN LONGARMI 
LONGARMI , LONGRMI / ADJCYC / PHIAI , PHI I / ADJCYC $ 

MERGBUS $ 

CHKPGEN $ 

QUIV3 ,NUPGEN $ 

CHKAGEN $ 

QUIV3 $ 

LONGPRMI , LONGRMI / NUPGEN / PH I P I , PH I I / NUPGEN $ 

MERGBUS $ 

CHKAGEN $ 

QUIV4 ,NUAGEN $ 

MERGBUS $ 

QUIV4 $ 

LONGARMI , LONGRMI / NUAGEN / PHIAI , PHI I / NUAGEN $ 

MERGBUS $ 

/, , , , 1-112 $ 

/LONGNULL, , , ,/0/2 $READ IN LONGNULL 
/LONGONE, , , ,/2/2 $ SKIP 2 DB & READ IN LONGONE 
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MERGE HMYPART, , , , ,L0NGNULL/LHRMHED/+7/+l $ INCREMENT = LENGTH OF NEW HARMONIC 
PARTN LONGONE, ,LHRMHED/ DUMMY, , ,/+7/+l $ PARTITION LONGONE DOWN BY HARMONIC 
MERGE DUMMY, , , , ,LONGNULL/HEDSHRTH/+7/+l $ APPEND HARMONIC-SIZE NULL TO TAIL 
ADD LONGONE, HEDSHRTH/LHRMTAL/ /( -1.0, 0.0) $HARMONIC-SIZE CLUSTER @ TAIL 

PARTN LONGRMI , ,LHRMHED/ DUMMY, , , / +7/+1 $TRIM A HARMONIC INCREMENT OF ZEROES 

$$ FROM THE HEAD OF LONGRMI. 

MERGE DUMMY, , , , , LHRMTAL /BBLRM/ +7 /+1 $PLUG A HARMONIC INCREMENT OF ZEROES 

$$ ONTO THE TAIL OF LONGRMI. 

ADD BBLRM , LHRMTAL / LONGRM J / $ CLUSTER =NUHARM INC + ACCEPTED 

PARTN LHRMTAL, ,LONGRMJ/ ,TRIMRM, , /+7/+1 $ THIS IS THE PARTITIONING VECTOR 
$$ NEEDED FOR COMBINING THE HARMONIC TO THE MATRIX OF ACCEPTED MODES 

$Test whether one or less columns of REMNANT are left. SET PARAMETER IF SO. 

$ 

PARAML REMNANT/ /*TRAILER*/l/V,N,REMCOL $ 

MERGE PHII , , REMNANT, ,TRIMRM, /PHIJ / +7/ +2 $ MERGED!! 

SWITCH LONGRMJ , LONGPRMI / /V,N,PRIMCYC $ 

SWITCH LONGRMJ, LONGPRMI/ /V,N,NUPGEN $ 

SWITCH LONGRM J,LONGARMI/ /V,N,ADJCYC 3 

SWITCH LONGRMJ , LONGARMI / / V , N , NUAGEN $ 

SWITCH PHIJ ,PHIPI/ /V,N,PRIMCYC $ 

SWITCH PHIJ ,PHIPI/ / V , N , NUPGEN $ 

SWITCH PHIJ ,FHIAI / /V,N,ADJCYC $ 

SWITCH PHIJ ,PHIAI / /V,N, NUAGEN $ 

SWITCH HEDVEC I , PHHED / /V,N,PRIMCYC $ 

SWITCH HEDVEC I, AHHED/ /V,N,ADJCYC $ 

SWITCH HEDVEC I, PHHED/ /V,N, NUPGEN $ 

SWITCH HEDVEC I, AHHED/ /V,N, NUAGEN $ 

SWITCH TAL VEC I , PHTAL / /V,N,PRIMCYC $ 

SWITCH TAL VEC I , AHTAL / /V,N,ADJCYC $ 

SWITCH TAL VEC I, PHTAL/ /V,N, NUPGEN $ 

SWITCH TAL VEC I, AHTAL/ /V,N, NUAGEN $ 

$ 
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PARAM / / *EQ* / V , N , HARMDONE/ C , Y , MODSPEC / V , N , HARMNO $ IF # HARM = MODS PEC =>DONE 

COND ORTHOG, HARMDONE $JUMF OUTSIDE HARMONIC LOOP IF HARMONICS ARE DONE 
$ 

COND PRIMOUT , PRIMCYC $ 

JUMP ADJHRMNY $ 

LABEL PRIMOUT $ 

SWITCH REMNANT, PGENI/ /V,N, PRIMCYC $ 

PARAM / /*MPY*/V, N, PRIMCYC /V,N, PRIMCYC/ -1 $ RESET PRIMCYC TO POSITIVE 

PARAM //*MPY*/V,N,ADJCYC/+1/-1 $ ENABLE THE LOOP FOR THE ADJOINT 1ST HARM 
JUMP HLOOPEND $ 

LABEL ADJHRMNY $ 

COND ADJOUT,ADJCYC $ 

JUMP PHMN2 $ 

LABEL ADJOUT $ 

SWITCH REMNANT, AGEN I / /V,N,ADJCYC 3 

PARAM / / *MPY* / V , N , ADJCYC / V , N , ADJCYC / -I $ RESET ADJCYC TO POSITIVE 

PARAM //*MPY*/V,N.NUPGEN/+i/-l $ ENABLE THE LOOP FOR THE PRIM HIGHER HARM 
JUMP HLOOPEND $ 

LABEL PHMNZ $ 

COND PHMNOUT , NUPGEN $ 

JUMP AHMNZ $ 

LABEL PHMNOUT $ 

SWITCH REMNANT , PGENI / /V,N, NUPGEN $ 

PARAM / /*MPY*/V,N,NUPGEN/V,N, NUPGEN/ -I $ RESET NUPGEN TO POSITIVE 

PARAM //*MPY*/V,N,NUAGEN/+1/-1 $ ENABLE THE LOOP FOR THE ADJ HIGHER HARM 
JUMP HLOOPEND $ 

LABEL AHMNZ $ 

SWITCH REMNANT, AGEN I / /V,N,NUAGEN $ 

PARAM //*MPY*/V,N,NUAGEN/V,N,NUAGEN/ -1 $ RESET NUAGEN TO POSITIVE 
PARAM //*MPY*/V,N,NUPGEN/+1/-1 $ ENABLE THE LOOP FOR THE PRIM HIGHER HARM 
LABEL HLOOPEND $ 

■5 


i 
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PARAM / /*LE*/V,N,TESTOVER/V,N,REMCOL/l $ TESTOVER =-l IF REMCOL <</ = ) 1 
$$ THIS IS THE EXIT IN CASE FEWER HARMONICS PASS 

$$ THE LI TEST THAN ARE SPECIFIED BY MODSPEC 

COND ALTCHK, TESTOVER $ EXIT WHEN HARMONICS ARE EXHAUSTED PRIOR TO MODSPEC 
JUMP USUAL $ 

LABEL ALTCHK $ 

COND USUAL, ADJCYC $ ADJOINT GETS A CHANCE TO GENERATE A SINGLE 
COND ORTHOG , NUPGEN $ PREVENT ANOTHER HARMONIC TO BE GENERATED FROM A SINGLE 
COND USUAL, NUAGEN $ ADJHRM GETS A CHANCE TO GENERATE A SINGLE 

JUMP ORTHOG $ 


LABEL 

PURGE 

REPT 

$ 

LABEL 

$ 

MERGE 

PARTN 

MERGE 

ADD 

PARTN 

MERGE 

PARTN 

PARTN 

LABEL 

OUTPUT1 , 

OUTPUT1 

LABEL 

PRTPARM 

LABEL 

END $ 


USUAL $ 

KMMTX/ MODSPEC $ 
HMNICGEN, 999 $ 


END OF HARMONIC GENERATOR 


r OOn o> o*, a, cv o. o, o^, o_ o„ o> ir, r>, r>, a, o, o, o, a, o_ c 


ORTHOG $ 


PHHED , , , , ,L0NGNULL/L0NGHEDl/+7/-H $**START OF HEADPHI CONSTRUCTION 

LONGONE, ,LONGHEDl/ DUMMY, , , /+7/+I $ LUMP TO MERGE ON HEAD 

DUMMY, , , , , LONGNULL/MISSTAIL/ +7 / +1 $ ONE MISSING FROM TAIL 

LONGONE, MISSTAIL/L0NGTAL1/ /( -1.0, 0.0) $**START OF TAILPHI CONSTRUCTIOI 

LONGONE, .LONGPRMI/ , DUMMY. ,/ + 7/+l $ALL ONES OF LGTH = ACCEPTED VECTORS 

DUMMY, , , , , LONGNULL /HEADER/ +7/+1 $HEAD CLUSTER OBVERSE OF LONGPRMI 

L0NGHED1 , .HEADER/ .HEADPHI, , /+7/+1 $ SAVE FOR DELIVERY TO ORTHOG 

LONGTAL 1 , .LONGPRMI/ .TAILPHI, , /+7/+1 $ SAVE FOR DELIVERY TO ORTHOG 

NOBIZNEZ $ GET OUT WITH OUTPUT SAME AS INPUT 

,,,,/ / — 1 / 3 $ SET DEFAULT LABEL 

PHIP I, HEADPHI, TAILPHI, PHI AI, / /0/3 $ 

PRINTOUT $ 

/ / 0 / C , N , HARMNO $ 

FINIS $ 

FINISH OF DMAP PROGRAM FOR RITZ HARMONICS 
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APPENDIX C 


RITZRTHG . DMP 

SELF AND DUO ORTHOGONAL I Z AT I ON 


NASTRAN MAXFILES = 60, FILES = (INP3,INP4) 

APP DMAP $ PROGRAMMED FOR 1988 NASTRAN. OUTPUT TO PUNCH FILE. 

$$ EXECUTES AFTER BOTH RITZFUND AND RITZHARM TO ORTHOGONALIZE RITZ MODES 
BEGIN $ $$ORTHOG. DMP 
INPUTT1 /, , , , /-1/3 $ 

INPUTT1 /PHIPI .HEADFHI ,TAILPHI ,PHIAI , /0/3 $ 

PARTN PHIPI, HEADPHI, / , ,PHI1, Z+7/+2 $ 

COPY PHI1/ZETA1/ 0 $ 

PARTN HEADPHI, , TA I LPH I / DUMMY , , , / + 7/+1 $ 

MERGE DUMMY HEADPHI /BBLHI/ +7 /+1 $ 

PARTN PHIPI, BBLHI, / , ,PHI2, Z4-7/+2 $ 

MPYAD PHI1 ,PHI2 , /NUM/+1/-1 $. 

COPY PHI 1 / CLONPHI 1 / 0 $ 

MPYAD CL0NPHI1 ,PHI1 , .'DEN/ + 1/+1 $ 

SCALAR NUM/ /l/I/ / /V,N,SPXNUM $ 

SCALAR DEN/ / 1/ 1/ / /V,N,SPXDEN $ 

PAR AMR / /*DIVC*/ / / /V.N,A11/V,N.SPXNUM/V,N,3PXDEN $ 

ADD PHI1 , FHI2/ZETA2/V,N,A11 $ SINGLE PREC. WON'T TAKE DEL PREC!!! 

PARTN TAILPHI, ,HEADPHI/DUMMD, , ,/ + //+-l $ 

MERGE DUMMD, , , , ,TAILPHI/BBLTI/+7/+l $ 

ADD BBLTI, TAILPHI /PTALCLUI/ $ 

PARTN TAILPHI, , PTALCLUI / .BUILDI , , /+7/+1 $ 

MERGE ZETA1 , ,ZETA2, , BUILDI, /ZETAI/+7/+2 $ 

MPYAD ZETAI, PHIPI, /C0EFI/+1/+I $ 

ADD BBLHI, HEADPHI /PHEDCLUI/ $ 

PARAM / /*ADD*/V,N,R0WC0W/2/0 $ 

PAR AML PHIPI/ /*TRAILER*/1/V,N,PC0L $ 
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LABEL ORTHLUP $ TOP OF SELF ORTHOGONAL I ZAT I ON LOOP 
PARAM / / *ADD* / V , N , ROWCOW/ V , N , ROWCOW / 1 $ 

PARTN COEFI , PHEDCLUI , / , ,CAI, /+7/+2 $. 

PARTN BBLHI , ,TAILPHI/DUMVEC , , , /+7/+1 $ 

MERGE DUMVEC, , , , , HEADPH I / BBLH J / +7/+1 $ 

PARTN COEFI, BBLH J, / , ,CFI, /+7/+2 $ 

SOLVE CAI , CFI /AIN/-1/-1/2 $ 

PARTN BBLTI , , HEADPH I / DM Y , , , /+7/+1 $ 

MERGE DMY, , , , ,TAILPHI /BBLTJ / +7/+1 $ 

ADD PTALCLU I , BBLTJ / PTALCLU J / $ 

PARTN TAILPHI , ,PTALCLUJ/ .BUILDJ, , /+7/+1 $ 

PARTN PTALCLUJ, .TAILPHI/ .UNIT. . /+7/+1/ 12 $ UNIT IS RECTANGULAR S.P. 

MERGE AIN, UNIT, , , ,BUILDJ/AJN/+7/+2 $ 

ADD PHEDCLUI . BBLHJ / PHEDCLUJ / $ 

PARTN PHIPI, PHEDCLUJ, / , ,PHIZ, Z+7/+2 $ 

MP YAD PHI Z , A JN , / ZETAX / 0 $ 

MERGE ZETAI , , ZETAX, , BUILDJ, /ZETAJ/+7/+2 S 

MPYAD ZETAJ, PHIPI, /COEFJ/+1/+1 $ 

SWITCH ZETAI, ZETAJ/ / -1 $ 

SWITCH BBLHI, BBLHJ/ / -1 $ 

SWITCH PHEDCLUI, PHEDCLUJ/ / -I $ 

SWITCH BBLTI, BBLTJ/ / -1 $ 

SWITCH PTALCLU I, PTALCLUJ/ / -1 $ 

SWITCH BUILDI, BUILDJ/ / -I $ 

SWITCH COEFI, COEFJ/ / -1 $ 

PARAM / / *EQ*/V,N,SELFDUN/V,N, ROWCOW/ V.N.PCOL $ 

COND DUALORTH , SELFDUN $ 

REPT ORTHLUP, 9 99 $ 

LABEL DUALORTH $ 

COPY ZETAI /CLONZETA/ 0 $ 

MPYAD CLONZETA, ZETAI, /ZSQ/+1/+1 $ 

MATPRN Z3Q ,,,,// $ 
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$ 

$ START OF DUAL ORHTOGONALIZATION OF ADJOINT MODES. 

$ 

INPUTT1 /, , , ,/-l/4 

INPUTT1 /MLL ,KLL , , , / 0 / 4 $ 

MPYAD ZETAI ,MLL, /ZEM/+1 $ 

MPYAD ZEM,PHIAI , /KOEF/O $ 

DIAGONAL MLL/UNITY/*SQUARE*/0.0 $ 

SOLVE KOEF, UNITY/ BETA/ -1/+1/+2/ +2 $ 

MPYAD PH I A I, BETA, / OMEGA /0 $ 

MPYAD OMEGA, MLL, /MEGM/+1/+1 $ 

MPYAD MEGM, ZETAI, /GENMASS/O $ 

MPYAD OMEGA, KLL, /MEGK/+1/+1 $ 

MPYAD MEGK, ZETAI, /GENSTIF/0 $ 

MATPRN GENMASS , ZETAI , OMEGA , GENSTIF , / / $ 

0 UTPUT3 ZETAI , OMEGA , GENMAS S , GENSTIF ,//0/C.Y,Nl=ZZZ/ 

C, Y,N2=MEG/C, Y,N3=MMM/C, Y,N4=KKK $ 

END $ FINISH OF ORTHOGONAL I Z AT ION OF RITZ VECTORS 
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APPENDIX D 


DEMONSTRATION PROBLEM FOR UNSYMMETRIC RITZ 
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APPENDIX E 


GENMASS KAPPA - .007 


1 . OOOOOE+OO 
7 . 33943E-25 
4 . 36876E-22 
6 . 08346E-27 


-1.57700E-21 
1. 00000E+00 
-1. 76725E-24 
-4.10067E-23 


-1 . 75392E-21 
5 . 52286E-26 
1. 00000E+00 
6 . 50197E-24 


-4.237 73E-26 
1 . 84410E-23 
-1 . 41962E-23 
1. 00000E+00 


GENMASS KAPPA =0.95 

1. 00000E+00 -3 . 41619E-17 -2.08085E-17 -7.77767E-17 -3.55682E-15 
6 . 80886E-17 7.19268E-16 -3.59782E-16 1.85463E-15 

-2. 09800E-17 1. OOOOOE+OO -2.66490E-16 5.34076E-15 4.56150E-17 

-1.0807 2732E-14 1.36521E-14 -1.98969E-13 

-7 . 27244E-16 3.14776E-15 1. OOOOOE+OO 1.58474E-15 6.56536E-16 
-7 . 71337E- 15 -3.17466E-14 1.96339E-14 -4.29536E-14 

-6 . 40719E-18 -1 . 66050E-15 5.80912E-17 1. OOOOOE+OO 4.28422E-19 

1 . 45089E-15 7 . 28284E-15 -5.87202E-15 1.86094E-14 

-1 . 15957E-15 1 . 05575E-16 4.21965E-17 1.36209E-16 1. OOOOOE+OO 
-2 . 06182E-16 -1 . 70354E-15 4.34321E-16 -2.34621E-15 

-4 . 61736E-17 4.20176E-15 -3.42951E-16 1.05352E-14 5.74604E-17 

1. OOOOOE+OO 2 . 34472E-14 1.11612E-14 -1.95622E-13 

-1 . 94247E-18 -4.76531E-16 3.61424E-17 -7.44025E-16 4.46676E-18 

1 . 22754E-15 1. OOOOOE+OO -1.59921E-15 2.45933E-14 

2 . 42787E-18 8.22934E-16 -2.48630E-17 1.54808E-15 -1.46353E-17 

-6 . 75527E-16 -7.19174E-16 1. OOOOOE+OO -1.18324E-14 

1 . 14403E-18 -1 . 60473E-17 -1.99198E-18 6.85738E-18 -1.30097E-18 

3 . 33071E-17 -6 . 84577E-18 -9.50914E-17 1. OOOOOE+OO 
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GENSTIF KAPPA = .007 


1 . 04682E+01 
-3 . 08635E-04 
-8 . 90200E-04 
-8 . 76229E-10 


6 . 15654E-04 
4 . 61758E-02 
-5 . 23407E-08 
2 . 54442E-05 


-9.26258E-04 
2 . 74302E-08 
1. 36253E-02 
-1 . 91582E-05 


-2.98338E-10 
5.87860E-08 
-3 . 94594E-08 
2 . 57602E+01 


GENSTIF KAPPA = 0.95 

3 . 39196E+02 -6.69263E-01 -2.44002E-01 -9.04235E-01 -3 . 12387E+02 
1.17373E+00 2 . 61338E+-01 -1.37057Ef00 2.05697E+01 

-1.08090E-02 3 . 15758E-01 -5.48427E-03 1.31833E-01 1.07321E-02 

-1.70708E-0I ‘ 8 . 53201E-01 6.03132E-01 -2.68652E+00 

-3 . 09520E-02 1.61339E-02 2.06822E-02 2.84852E-02 2.85713E-02 
-3 . 69205E-02 -2.37364E-01 1.00253E-01 -6.03664E-01 

-8 . 47928E-01 1.11503E+02 -8.23610E+00 6.49374E+02 1.19001E+00 

-8 . 02691E+02 3.80773E+02 3.22937E+02 -4.57428E+03 

-3 . 07444E+02 6.23311E-01 2.25336E-01 8.38877E-01 2.89818E+02 
-1 . 08889E+00 -2 . 42260E+01 1.27283E+00 -1.90819E+01 

7 . 98571E-01 -1 . 04804E+02 7.74130E+00 -5.92276E+02 -1 . 12001E+00 

7 . 54470E+02 -3.57896E+02 -3 . 03535E+02 4.29948E+03 

1 . 98390E-01 2 . 15040E-02 -8.83532E-03 2.23315E-02 -1.85797E-01 

-2.80937E-02 3.86155E-01 -1.58374E-02 6.00169E-02 

6 . 36402E-02 1.95278E+00 -1.21436E-01 4. 85570E+00 -5.51377E-02 

-6 . 20209E+00 6.10810E+00 5.60222E+00 -7.46195E+01 

1 . 76325E-02 -2.24893E+00 1.66211E-01 -5.93498E+00 -2.44872E-02 

7 . 57854E+00 -7.62053E+00 -6.48352E+00 9.22406E+01 
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